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Abstract 

We apply phase-space density considerations to obtain lower bounds on the mass of 
sterile neutrino as dark matter candidate. The bounds are different for non-resonant 
production, resonant production in the presence of lepton asymmetry and production 
in decays of heavier particles. In the former case our bound is comparable to, but 
independent of the Lyman-a bound, and together with X-ray upper limit it disfavors 
non-resonantly produced sterile neutrino dark matter. An interesting feature of the 
latter case is that warm dark matter may be composed of heavy particles. 

One candidate for dark matter particle is sterile neutrino. In the first place, the existence 
of sterile neutrinos is favored by the discovery of neutrino oscillations. Furthermore, by 
adding sterile neutrinos to the Standard Model and fine-tuning parameters in the neutrino 
sector, one is able to explain all known facts in high energy physics and cosmology without 
introducing any other new fields except for inflaton [1, 2, 3]. 

The sterile neutrino dark matter should satisfy cosmological and astrophysical con- 
straints. These constraints depend on the way sterile neutrinos are produced in the early Uni- 
verse. There are several mechanisms capable of generating sterile neutrino in right amount. 
These include the non- resonant production via active-sterile neutrino mixing [4], resonant 
production in the presence of lepton asymmetry [5, 6], production in decays of heavier par- 
ticles (see, e.g., Refs. [3, 7] in the context of z/MSM model) and production in scattering 
(see e.g., Ref. [8]). In this paper we consider astrophysical constraints on sterile neutrino 
dark matter generated by the non-resonant production mechanism. We comment on other 
production mechanisms towards the end of this paper. 

Sterile neutrino is coupled to active ones through mixing only. Thus, the mixing angle 
controls the decay rate of sterile neutrino into active neutrino and photons. These decays 
in galactic halos would produce X-ray emission lines, which are severely constrained [9] by 
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current observations. Within the non-resonant production mechanism, the mixing angle 
and the mass m of sterile neutrino are related to each other by the requirement that sterile 
neutrinos make all of the dark matter. So, the X-ray observations put the upper bound on 
the mass, m < 4keV [10]. 

Sterile neutrino with the mass in the keV range is warm dark matter candidate. In fact, 
warm dark matter may be even more attractive than cold dark matter [11], as it might 
be able to address emerging problems of the standard CDM cosmology on small scales, 
such as missing satellites [12], galactic density profiles [13] and angular momentum of spiral 
galaxies [14]. However, like any warm dark matter candidate sterile neutrino should not be 
too "warm" in order to form observed small scale structures. This implies that its mass 
should not be too small. 

Well known lower limits on sterile neutrino mass come from the observations of Lyman-a 
forest — multiple absorption lines in spectra of distant quasars [15, 16]. Most of these lines 
correspond to smooth overdense regions of warm ionized intergalactic medium which are 
believed to trace dark matter clustering. To derive the sterile neutrino mass bounds, one 
compares the statistical properties of Lyman-a absorption lines with those predicted within 
various cosmological models by semi-analytical methods and numerical simulations. Recent 
Lyman-a lower bounds on the mass of sterile neutrino dark matter produced via the non- 
resonant production mechanism are m > 5.6 — 28keV depending on data set [16]. Together 
with the X-ray upper bound they disfavor the non-resonantly produced sterile neutrino as 
the dark matter candidate. 

Given the uncertainties in the Lyman-a constraints, it is worth deriving a lower bound on 
sterile neutrino mass in an independent way. In this paper we employ the Tremaine-Gunn 
approach, which was successfully used nearly thirty years ago to rule out the scenario of 
active neutrino hot dark matter [17]. This approach relies on the fact that coarse-grained 
distribution function of collisionless particles decreases in the course of evolution. 

The process of gravitational formation of compact objects from collisionless dark mat- 
ter particles is described by the Vlasov equations governing the evolution of distribution 
function in self-consistent gravitation field, see, e.g. Ref. [18] and references therein. One 
of the most important phenomena taking place in such systems is mixing of distribution 
function in phase-space. The phase-space volume occupied by the system is invariant under 
Vlasov equation dynamics. Nevertheless, particle trajectories chaotically disperse, explor- 
ing new areas of phase-space and forming there ever more complicated fine structure. This 
results in overall decrease of coarse-grained distribution function. In particular, values of 
coarse-grained distribution function cannot exceed the maximum value of the primordial 
distribution function [19]. The latter property enables one to obtain constraints on dark 
matter. 

Primordial distribution function of sterile neutrinos produced via non-resonant oscilla- 
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tions can be approximated by scaled distribution of active neutrinos — ultra-relativistic 
Fermi-Dirac distribution [4], 

rip) = -1 L_ (i) 

J[P) (2tt) 3 eP/^ + 1 ' K > 

where g = 2 is the number of sterile neutrino spin states. The coefficient (5 here is propor- 
tional to the active-sterile mixing angle squared and can be related to the sterile neutrino 
mass by demanding that sterile neutrinos constitute all dark matter in the Universe 1 : 

'^dm\ /10eV s 



1 0.2 J V m 

Thus, the maximum value of the primordial distribution function in this scenario is given by 

max/(p ) 1 (2) 
J W (2tt) 3 V 0.2 J V m J K J 

Observationally, the value of the coarse-grained distribution function in a galactic halo 
can be estimated by the phase-space density, the ratio between the mass density and cube 
of the one-dimensional velocity dispersion in a given volume, Q = p/a 3 [21]. In practice, 
one measures the velocity dispersion of stars and assumes that it coincides with that of dark 
matter particles. For non-relativistic dark matter particles one has 

n 



Q = m 



i 



(l p 2}3/2 ' 



where n is their average number density in a halo. Assuming that the coarse-grained distri- 
bution of halo particles is isotropic, fhaio(p, r ) = fhaio(p,r), one estimates 



n 



5/2 



[I halo (P,r)rf 3 p] 

3/2 
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fhaioip*,r) , (3) 



where is a typical momentum of the dark matter particles. 2 In this way the value of the 
coarse-grained distribution function in a galactic halo is estimated as 

f hal ° ~ 33/2 m 4 • ( 4 ) 
Hence, one arrives at the following constraint on dark matter, 

Q 



33/2 m 4 



< max f(p) . (5) 



1 We do not consider here the possibility that the non-resonant mechanism generates not all, but only a 
fraction of dark matter [20]. 

2 In the case when the width of the momentum distribution around p* is small, Ap < p*, the estimate (3) 
reads n/(p 2 ) 3 / 2 <~ fhaio • (Ap/p*). Then, instead of (4), one has an inequality fhaio > Q/(3 3 / 2 to 4 ). This 
makes the constraint (5) even stronger. 
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The strongest constraints on warm dark matter scenario are obtained by making use 
of the highest values of the phase-space density in dark matter dominated objects. These 
are found in dwarf spheroidal satellite galaxies [22, 23]. dSph's are the most dark matter 
dominated compact objects observed so far, and are conjectured to be hosted by the smallest 
possible dark matter halos [22]. In recently discovered objects Coma Berenices, Leo IV and 
Canes Venaciti II, the value of Q ranges from 5 • 1(T 3 to 2 • 10~ 2 [24]. We use 

' ^ to (km/s) J (km/s)' 1 L 1 

the first, slightly more conservative value, 

0s9 ^ = 5.i(r«f^. (6) 

(km/s) (km/s) 
Making use of Eqs. (2) and (5) we obtain the constraint on the sterile neutrino mass, 

->"-0' /3 (a^) 1/3 - 

Despite the fact that the value of observable Q may not be exactly the same as the 
value of the coarse-grained distribution function of dark matter particles, we consider the 
bound (7) as conservative. Indeed, the coarse grained distribution function is likely to de- 
crease considerably during the non-linear stage of evolution. For example, some numerical 
simulations of halo formation show the decrease of Q by a factor of 10 2 — 10 3 from input 
initial values [25]. Therefore, further improvement of understanding of how compact objects 
are formed by warm dark matter particles is likely to strengthen the bound (7). We con- 
clude that the Tremaine-Gunn approach gives another argument, independent of Lyman-a 
, that disfavors the sterile neutrino dark matter generated by the non-resonant production 
mechanism. 

However, sterile neutrino may still be a dark matter candidate. There are other generation 
mechanisms of sterile neutrino dark matter, such as resonant production in the presence of 
lepton asymmetry [5, 6], production in scattering [8] and production in decays [3, 7]. In 
these cases the sterile neutrino decay rate is not directly related to their abundance in the 
Universe, so the X-ray mass bound does not apply. 

Although the primordial distribution function of sterile neutrinos in these scenarios differs 
from (1), lower bounds from both Lyman-a and phase-space observations are not expected 
to change dramatically. For example, the results presented in Ref. [6] indicate that the 
maximum value of the distribution function of 3 keV sterile neutrinos generated by the 
resonant production in the presence of lepton asymmetry is about 0.03 — 0.3 of the maximum 
of thermal distribution, l/(27r) 3 (accounting for 2 spin states), depending on the value of the 
lepton asymmetry. In that case one has 

3 3 / 2 m 4 max/Go) = (6.5 • lO" 3 - 6.5 • lO" 2 ) , 

(km/s) 
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at the edge of the constraint (5), with Q estimated as in (6). In this case the correct dark 
matter density, Q u ~ 0.2, is obtained for the mixing angle in the range sin 2 (26) = 10~ n — 10~ 8 
and even smaller, again depending on the lepton asymmetry, and the X-ray bound can be 
satisfied [6]. 

If sterile neutrinos are produced in scattering processes and do not equilibrate, as is the 
case for the production mechanism of Ref. [8], then the neutrino distribution is approximately 
given by (1) with the same coefficient /3 and effective temperature T v . Hence, the maximum 
of the neutrino distribution function is the same, as in the case of non-resonant production. 
So, the limit (7) is applicable for sterile neutrino produced in scattering as well. 

Before proceeding to sterile neutrino production in decays, let us recall that there is a 
simple bound on their mass based on Pauli blocking. Namely, the primordial distribution 
function cannot exceed the value 2/(27r) 3 , so for distributions saturating this bound, the 
constraint (5) translates into 

m > 1.0 keV ( g ) V4 , (8) 



5-10" 

where we again used the estimate (6). This is a model-independent lower limit on the mass 
of a fermionic dark matter candidate (assuming two spin states). In what follows we do not 
take into account Pauli blocking, with understanding that if the limits obtained are weaker 
than (8), then the limit (8) applies. 

Sterile neutrinos may be produced in decays of relativistic thermalized particles 3 of mass 
M and partial decay width at rest into sterile neutrinos T. In that case the low momentum 
part of sterile neutrino distribution function is given by [26, 3, 7, 27] 



3 M 2 \ p J J v ' 4tt 5 /2 M 2 \ p 

o 

Here Mp X = M m y/90/(87v 3 g*), T 0jeff = (^f^J T , where g* and g* fi are the effective 
numbers of degrees of freedom at decay and present epoch, respectively, and f t h is the 
thermal distribution function of decaying particles. We assume here that the latter are 
scalars. Hereafter f(jp) denotes the primordial distribution function redshifted to the present 
epoch. The total present number density is [26] 

_3C(5) 3 M^r 



3 Wc assume here for simplicity that these particles interact sufficiently weakly with the rest of the cosmic 
plasma, so that sterile neutrino production in scattering processes is negligible. In the opposite case the 
mass bound is similar to (7). 
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Requiring that sterile neutrinos make all of dark matter, n^m = QdmPc one finds the only 
relevant combination of parameters of decaying particles, 



r 47T ^DM Pc 

M2"3C(5) mM* x Tl ff - 

The distribution (9) is formally unbounded from above. In reality this means that at low 
momenta, the distribution function takes the Pauli blocking value, / = 2/(27r) 3 . In this 
situation one in principle is still able to obtain mass limits stronger than (8) by invoking 
the following statistical argument [28]. One requires that in the early Universe, a certain 
fraction v of dark matter particles are sufficiently densely packed in phase space so that 
these particles are able to form subsequently dark matter halos of high Q. In other words, 
the value of the distribution function of this fraction of particles should obey the constraint 

m > 33^4 , no) 

where the observed value of Q is estimated as in (6). The fraction v should not be smaller 
than the fraction of dark matter residing in dSph's, which we estimated in [26] as 

v ~ KT 5 . 

Given very weak dependence on v in the limits we obtain in what follows, the precise number 
is unimportant for our purposes. 

We continue the discussion of thermal creation by noticing that the fraction v of most 
densely packed particles is related to the maximum momentum p v of these particles in an 
obvious way, 

f(p) Anp 2 dp = un . 



o 



Making use of (9) we obtain 



Pf 



i5v^c(5) y /5 



T , e// V 8C(5/2) 
Substituting this value of momentum back into (9) we find 

; ^-mT 3 e// 3^U(5) J U; 



1.1 • 10 -2 



^dm\ / 9* \ /lkeV\ /10" 5 V /5 



0.2 J V106.75/ V m J \ v 
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It is this value that should obey the constraint (10). Hence, we obtain the lower bound on 
the sterile neutrino mass, 

m > 0.88keV f^V' 3 f^V /3 f^m 1 ' 3 (A)" 15 . (11) 
V^dm/ V 9* J V5-10- 3 / VlO" 5 / V ; 

This bound is in fact slightly weaker than the Pauli blocking bound (8) for our values of 
parameters. As we noticed above, in this situation one should use the bound (8) instead. 

The bound (11) is somewhat different from the bound obtained in Ref. [27], since we use 
the statistical approach rather than the approach of Ref. [21]. Notice that the dark matter 
is warm, i.e., the bound (8) is nearly saturated if the parameters M and T are such that 

M^~10-GeV /02W106 75X3 ' 2 



r y Odm /V 9* 

This is the case if either the decaying particles are very heavy or their decay rate into sterile 
neutrinos is very small, or both. 

Finally, let us consider the production of sterile neutrinos in decays of non-relativistic 
particles whose number in comoving volume has been frozen out. In that case the momentum 
of sterile neutrinos at production equals = M/2, and then the momentum gets redshifted 
to 

P = P* -JTT = P*H^ , 12 

a(t ) T 

where t is the time at decay and to is the present time. Hence, the distribution function of 
sterile neutrinos is obtained from the relation 

f(p) d i P = n e- v ^ t T tot dt. 

Notice that we normalized the distribution function to the present number density of sterile 
neutrinos, assuming that the abundance of decaying particles is just right to produce them. 
Therefore, this formula contains the total width Y tot only. Making use of (12) we obtain 

/(pHnoe -r,^J_. (13) 

Now, at radiation dominated epoch 4 t = Mpj/(2T 2 ), so that from (12) we obtain 

Mpj fpV 1 



^^uAt) - (14) 



4 One can check that the energy density of decaying non-relativistic particles never dominates in the case 
we consider here. 
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Hence, the distribution function (13) behaves as 1/p at low momenta, and we again have 
to employ the statistical argument. Time t v by which the fraction v of sterile neutrinos is 
produced, is determined by 

v — i — e — i tot l u 



The corresponding momentum p u is found from (14). Inserting it into (13) and requiring 
that sterile neutrinos make all of dark matter, we find 

y/2 I^dm Pc f M^Tto^ 3 ^ 2 _ 1/2 



rr nnl ff \ M 2 



2 4 f M Pi r tot \ 3/2 (n DM \ f g* \ /lkeV\ / 10^ \ 1/2 
' \ M 2 ) V °- 2 / V106.75/ V m J \ v ) 

We again make use of (10) for this value of the distribution function, and obtain finally the 
bound 

M 2 \ 1/2 ( 0.2 V /s ^106.75 V /3 ( q ^V3 ( v \V6 



m>145eV feJ 'fej {—) fe) ^) < 15 > 



This bound should be used whenever it supersedes the bound (8). 

Unlike the limit (11), the bound (15) depends on parameters M and T tot characterizing 
the decaying particles. This is because unlike in the thermal production case, we now have 
(implicitly) one more free parameter, the number density of decaying particles at their freeze- 
out. Interestingly, for heavy enough decaying particles and/or long enough lifetime of these 
particles, the right hand side of (15) may be well above the "canonical" keV range. This 
implies that the decay mechanism we discuss here is capable of producing warm dark matter 
composed of heavy particles. As an example, for M ~ 10 14 GeV, Y tot = f^M and y ~ 10~ 12 , 
sterile neutrino of mass in TeV range would be warm. 

It is worth noting that the bounds (11) and (15) apply not only to sterile neutrinos but 
to any fermionic dark matter candidates produced in similar decay processes. 

After this work has been completed, we received a draft of the paper [29], where similar 
issues have been considered. Our results are consistent with the results of Ref. [29] wherever 
they overlap. 
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